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ABSTRACT 

We report on two XMM-Newton observations of the planetary host star 
HD 189733. The system has a close in planet and it can potentially affect the coro- 
nal structure via interactions with the magneto sphere. We have obtained X-ray spectra 
and light curves from EPIC and RGS on board XMM-Newton which we have analyzed 
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and interpreted. We reduced X-ray data from primary transit and secondary eclipse oc- 
curred in April 17th 2007 and May 18th 2009, respectively. In the April 2007 obser- 
vation only variability due to weak flares is recognized. In 2009 HD 189733 exhibited 
a X-ray flux always larger than in the 2007 observation. The average flux in 2009 was 
higher than in 2007 observation by a factor of 45%. During the 2009 secondary eclipse 
we observed a softening of the X-ray spectrum significant at level of ~ 3o~. Further, 
we observed the most intense flare recorded at either epochs. This flare occurred 3 ks 
after the end of the eclipse. The flare decay shows several minor ignitions perhaps 
linked to the main event and hinting for secondary loops that emit triggered by the 
main loop. Magneto-Hydro-Dynamical (MHD) simulations show that the magnetic 
interaction between planet and star enhances the density and the magnetic field in a re- 
gion comprised between the planet and the star because of their relative orbital/rotation 
motion. X-ray observations and model predictions are globally found in agreement, 
despite the quite simple MHD model and the lack of precise estimate of parameters 
including the alignment and the intensity of stellar and planetary magnetic fields. Fu- 
ture observations should confirm or disprove this hypothesis, by determining whether 
flares are systematically recurring in the light curve at the same planetary phase. 

Subject headings: stars: activity - planetary systems - stars: magnetic fields - stars: 
individual (HD 189733) 



Introduction 



Recent work indicates that exoplanets, and especially Jupiter class planets at a distance < 0.1 
AU from their parent stars (the so called hot Jupiters), are in a unique X-ray environment and 
the effects of stellar X-rays may be significant to their evolution. The X -rays have been cite d 
as the cause of excess heating of the planet, which can induce mass loss (|Lammer et al.ll2003|) . 
Further, it has been argued that the magnetic fields of the star and the planet can interac t when 
the latter is orbiting at a distan ce of very few stellar radii as in the case of HD 189733 (|Lanza 
2008L |2009j; ICohen et al.ll2009|) . The complex magnetosphere that is formed by the star+planet 
system can drag material from the outer atmosphere of the planet and funnel it onto the stellar 
surface. The accretion of gas on the star that eventually could arise from this phenomenon can 
provoke shocks that heats the plasma up to a few million of degree and thus emits in X-rays. 
Furthermore, the enhanced magnetic field near the stellar surface can form very active regions 
on the star increasing the activity of the star in the X-ray band. Recent statistical analysis of X- 
ray activity of stars possessi ng hot- Jupiters ind icates that stars with close orbiting planets show 



enhanced X-ray emission (|Kashyap et al.ll2008|) . In the optical band. IShkolnik et al 



(2003) have 
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detected enhanced activity of Ca H & K in phase with planet's orbital period in HD 179949. 



A typical "hot Jupiter" receives ~ 2 x 10 4 times more radiation from its star than Jupiter 
does from the Sun, and orbits at a distance of only ~10 stellar radii. Unlike solar system planets, 
where residual heat from formation and gravitational settling still plays an important part in the 
energy budgets of Jupiter and other giant planets, the energy budget for hot Jupiters is completely 
dominated by the strong radiation that they receive from their parent stars. This energy budget 
is poorly constrained in the X-ray and UV regimes by standard stellar atmosphere models. Thus, 
understanding the nature of this input and the planet-star interacti ons is essential for un derstanding 
the properties and the evolution of these planets. As pointed by iLammer et al.l (|2003|) . the stellar 
X-ray and UV fluxes in their irradiance calculations strongly determine the amount of losses of 
gas from the outer atmospheres of close-in exoplanets, otherwise the losses of volatiles in models 
neglecting UV/ X-ray irradiation are proved to be unrealistically low. The UV/X-ray flux from 
the star plays a crucial role in determining the photochemistry of the planet's upper atmosphere. 
The photo-chemical products can act, in turn, as high-altitude absorbers, create thermal inversions, 
and otherwise alter the observable properties of the planet's upp er atmosphere. These effects are 



impo rtant when observing at lower energies/longer wavelengths (Bu rrows et al.ll2008l iLiang et al. 
20041) 



HD 189733 (R.A. 20 /? 00 m 43.7*, Dec. +22 rf 42 m 39.U) is a K 1.5V type star at a distance of 
193 pc, of mass 0.81 Mo, radius 0.76 R© and a rotational period of ~ 11. 95 d (|Henry & Winn 



2008J). It is in a wide binary system and HD 189733B, the oth er stellar member of the binary 
system, is a M4V type star (with range: M3.5-M5, M ~ 0.2M Q , Bakos et al.ll2006|) orbiting at a 
mean distance of ~ 220 AU (12.3" apparent separation) with a period of 3200 yr in a plane nearly 
perpendicular to the Earth-HD 1 89733 A line of s ight. It is thought that the age of this system is 
greater than or equal to 0.6 Gyr (IMelo et al]|2006|) . The magnetic field is estimated of a few tens 
of Gauss, a toroidal component and higher multipole components are presumed and a strong shear 
of field lines due to differential rotation is inferred (Mo utou et al.ll2007 ). 



HD 189733 A hosts a Jupiter size planet (HD 189733b, M ~ 1.15 M Jup , R ~ \ 26R iup ) at a 



distan ce of only 0.031 AU (8.75 stellar radii) and with an orbital period of 2.219d (|Bouchy et al 



2005J), much lower than the stellar rotational period. HD 1 89733b is thus classified as a hot Jupiter 
planet type. The orbit is nearly circular and orbital plane is parallel within a few degrees of the 
Earth-primary line of sight. This system has been the subject of many studies since its discov- 
ery, given its close distance and favorable geometry. Given Q-values for giant planets in the solar 
system, it is thought that orbital radii of a few hundredths of an AU from parent star reduces hot 
Jupiters tidal dissipation timescales to a few Myr, and tidally locked hot Jupiters should be the 
norm. Under this hypothesis, HD 189733b should show a pr onounced temp erature difference be- 



tween the day-side and the night-side hemispheres. Instead iKnutsonl (|2007|) found only a modest 
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difference of temperature between day-side and night-side hemispheres, suggesting that the ther- 
mal energy transport and thermalization in the planet's atmosphere is quite efficient. The same 
study shows that a bright spot is present in the planet surface displaced by 16 + 6° east longitude 
with respect to the sub-stellar point. 



Observations with ACS camera on board Hubble allowed lPont et al.l (|2007|) to infer that planet 
transits cover spots on the stellar surface with sizes of order of 80000 km. HD 1 89733 A is known to 
be an active and varia ble K type star with strong magnetic activity, as inferred from Call H&K lines 
(IMoutou et al.l 120071) . Through modeling the observed Rossiter-McL aughlin effect (the apparent 
Doppler shift occurring for the transit of the planet at the stellar limbs') IWinn et al.l (|2006) find that 
the spin stellar axis and the orbit normal are misaligned of only 1.4", hinting that migration of the 
planets does not change this alignment. 

Due to its close distance and the favorable orientation of the planetary orbit plane, HD 1 89733 
is an ideal laboratory for the study of close in planets. Among the plethora of known exoplanets, 
the details on the atmosphere of HD 189733b obtained through optical and IR spectroscopy are 
unique. 

In this paper we report on two observations of HD 1 89733 system obtained in X-ray band with 
XMM-Newton satellite at the time of a planetary transit and a secondary eclipse, respectively. X- 
ray observations can constrain model predictions, reduce degeneracy in their parameter space and 
help MHD simulations to infer physical and geometrical properties of the complex magnetosphere 
created by star and planet. 

The structure of the paper is as following: in Sect. [2] we describe the observations and the 
data analysis; In Sect. [3]and|4]we show the results and discuss them and in Sect. [5]we present our 
conclusions. 



2. Observations and Data Analysis 

We have obtained an observation of HD 189733 with XMM-Newton during the eclipse of 
the planet (secondary transit, when the planet passes behind the star as seen from the Earth) that 
occurred on 18 - 19 May 2009 (Obs # 0600970201). We have also analyzed a previous XMM- 
Newton archival observation of HD 189733 during the planetary transit of 17 April 2007 (Obs. 
# 0506070201, P.I. P. Wheatley). A source close to HD 189733 is visible in both observations, 
especially in 2007, while the M-type companion HD 189733B is undetected. Obviously, the planet 
is not resolved. Fig. [T] shows the EPIC images of HD 189733 obtained in 2007 (top row) and 2009 
(bottom row), respectively. The XMM Observation Data Files of both observations were processed 
with SAS (ver. 8.0) to obtain event tables of MOS1, MOS2 and PN calibrated both in energy and 
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Fig. 1 — EPIC images of HD 189733 A+B observed with XMM-Newton in 2007 April (top row) 
and 2009 May (bottom row). Images of MOS 1, MOS 2 and PN are shown in sequence from left 
to right. The companion HD 189733B is not visible in both observations. A bright source with 
uncertain optical counterpart is visible to the south of HD 189733A, and it has faded in 2009. 
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HD 189733A - Ml 



HD 189733A - Ml 




Time (ks) since 2007-04- 1 7T1 4: 1 1 : 1 9. 1 63 



10 15 20 25 

Time (ks) since 2009-05- 1 8T2 1 :41 :01 .093 



Fig. 2.— Light curves of HD 189733. On the left column: light curves of MOS1, MOS2 and 
PN acquired during the 2007 observation. On the right column the same for the 2009 observation. 
Error bars are at lcr level. The red light curve is in the 0.3-0.8 keV energy band while the blue 
light curve is in 0.8-2.5 keV energy band. Gray curve is the background. Bin size is 900 s for PN 
and 1800 s for MOS 1 and 2. 
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position. We restricted the analysis to events in the 0.3 - 8.0 keV energy band and to events that 
triggered at most two adjacent pixels as suggested by the SAS guide. We have extracted events 
of HD 189733A from a circular region that avoids the contamination from the M star companion 
HD 189733B. From the events selected in this region we have obtained light curves (Fig. [2]) and 
EPIC spectral of HD 189733A. For the analysis of light curves and spectra we focus primarily on 
the PN data in this work because of its higher efficiency in collecting counts with respect to both 
MOS cameras (4557 PN counts vs 1166 MOS1 and 1190 MOS2 counts in the source extraction 
region, respectively). A likely loss of data telemetry causes a dip at ~ 23 ks with a duration of 
~ 900 s in the light curves of 2009 observation, more marked in PN data. 

In order to characterize the variations of the source at short timescales, we consider the pulse 
invariant (PI) spectrum of PN events in a moving window comprised of 200 events. We compute 
the average energy of the spectral distribution, as well as the median, the standard deviation, and 
the 16% and 84% quantiles for a contiguous set of 200 photons (see Figure©. The window 
is moved by adding five events at the end and removing five from the beginning. We obtain a 
smoothed light curve in the same manner, with the time interval between the first and last events 
in the window suitably corrected for the presence of gaps in the exposure time. A background is 
estimated for each interval and is subtracted from the light curveQ. 

We have also extracted RGS spectra to have a higher resolution spectrum (SA ~ 0.04A) in 
order to observe the behavior of features, especially O VII triplet and its changes during the ob- 
servation. The ratio of forbidden to recombination lines from OVII triplet offers a diagnostic of 
the plasma density. For the main flare observed in the light curve of 2009 observation, we have 
performed time resolved spectroscopy of PN data by considering the spectrum during the flare, 
as well as before and after the flare. The spectra of EPIC PN have been analyzed with XSPEC 
software to obtain best fit parameters from thermal models (APEC) with one or two temperatures. 
Given that absorption is very low due to the short distance to the star, we have considered only 
unabsorbed spectral models. The free parameters were the temperatures, the global abundance and 
the emission measures. 



3. Results 

In Fig. [2]we plot the light curves from the 2007 and 2009 observations in three bands: 0.3-0.8 
keV (red), 0.3-2.5 keV (blue) and 0.3-8.0 keV (black). The background curve is shown in gray. 
The duration of the planetary transit and the secondary eclipse is ~ 108 min from 1st to 4th contact, 



'Note that background is ignored for the spectral distribution both due to the relative constancy of the background 
spectrum and the numerically small correction to the spectrum; see Figure|2] 
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Obs. 2009: TIME [ + 359068. 8,ks] 
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Fig. 3. — Intensity (top panels) and spectral variation (bottom panels) of HD 189733 A during 2009 
(secondary eclipse) registered from PN. On the left panels the curves relative to the whole exposure 
time. Vertical lines denote the 1st, 2nd, 3rd and 4th contact. In the bottom left panel, the green 
curves represent the 16% and 84% quantiles of the distribution of a sample of 200 photons, and 
show the overall width of the PI spectral distribution. Red curves represent the mean ±lcr standard 
deviation of the mean of the sample of PI between 0.3 and 1 .5 keV, and show the temporal variation 
of the spectrum. On the right panels we plot the rate and the mean of PI during the flare with a 
log scale in the y-axis of rate. Vertical lines denote the beginning of the flare, the peak, and three 
secondary impulses. 
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Viewing geometry at phase 0.5 



Viewing geometry at phase ~0.54 
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Fig. 4. — Scheme of the angle between sub-planetary point and active region where the large flare 
could have occurred in 2009. 
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Fig. 5. — Left panel: PN spectra of HD 189733 A during the three phases, pre-flare (0-15 ks, 
black), flare (15-25 ks, red) and post-flare (25-35 ks, blue). Right panel: flare spectrum (minus 
pre-flare spectrum) with best fit model and residuals. The region of OVII at ~ 0.57 keV (~ 22A) 
is enhanced during the flare. 
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RGS1 — pre flare 
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Fig. 6. — RGS 1 spectra of HD 189733A during the three phases: pre-flare (0-15 ks), flare (15- 
25 ks) and post-flare (25-35 ks). 
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i.e. 6480 s. 

The 2007 light curve shows a modest flare in coincidence of the planetary transit as the spec- 
trum seems to harden at ~30 ks (cp ~ 0) and decays exponentially post-transit. We find that a 
1-T thermal model with plasma temperature of ~ 0.6 keV reasonably fits the spectrum during the 
first flare. This is the most prominent feature in the 2007 light curve. The remaining part of the 
observation is best modeled with a 2-T thermal model with a cool component (kT = 0.24 keV) and 
a hot component (kT = 0.56 keV) very similar to the first flare temperature. The unabsorbed flux 
in the 0.3-8 keV band energy is ~ 2.5 x 10~ 13 erg s" 1 cm -2 for a luminosity of L x = 1.1 x 10 28 erg 
s" 1 . Table \T\ shows the results of the fits. Abundances ar e found at Z ~ 0.15 - 0.2Z G , as in active 



stellar coronae and very young stars (Maggi o et al.ll2007|) 



During the 2009 observation the star had a X-ray activity more pronounced than in 2007. The 
mean energy of the 2009 spectrum is ~ 700 keV and the PN rate is on average ~ 65 ct ks" 1 during 
most of the observation, reaching a peak of ~ 95 ct ks" 1 during the first flare. In 2009 the PN rate 
of HD 189733A is always larger than in 2007: the lowest rate is in fact ~ 95 ct ks" 1 in 2009 with 
a peak of ~ 350 ct ks" 1 (see Fig. [3]). 

During the secondary eclipse, the light curve shows a softening of the spectrum. From Fig. |3] 
we see that the mean of the PN distribution of PI gets softer during the eclipse. The mean photon 
energy in the 0.3 -1.5 keV range before and after the eclipse is ~ 700 + 10 (l<x), while during the 
middle of the eclipse it is ~ 660 + 10 eV (lcr), thus it is significant at a level of 3 ~ cr. During 
the eclipse the rate increases. The quartile curves show that the soft part of the spectrum increases 
while the hard part of the spectrum decreases. 

The light curves in Fig. |2] and |3] show a flare with a decay of ~ 10 ks that occurs 3 ks after 
the end of the planetary eclipse. The flare decay is characterized by several smaller impulses after 
the main ignition^. The flare phase is cp ~ 0.54 which would correspond to the planet being 16° 
past the earth-star line of sight (see Fig. |4]). In the soft band (0.3-0.8 keV) the light curve of PN 
has a step-like profile (Fig. [2]), characterized by an almost flat high level whereas in MOS 1 and 
2 the soft band light curves show a more usual peaked profile. The difference of the flare profile 
in PN and MOS is very likely due to their different responses at low energies with the PN being 
more sensitive to low energies than MOsO The plasma spectrum becomes hard at the flare peak 
as shown by the median and quartile curves (right panel in Fig. [3]). We deduce that the temperature 
of the plasma increases approaching the peak 



2 A likely loss of telemetry is causing the sharp dip in the decay around 23 ks. 

3 The ratio of normalized PN and MOS effective areas shows that PN is proportionally more sensitive than MOS 
by a factor ~ 2 in the range 0.3 - 0.7 keV. 
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To investigate the characteristics of the flare we have analyzed the PN spectra in the intervals 
0-15 ks, 15 - 25 ks, and 25 - 35 ks (chosen by visual inspection of the light curve). In TableCQwe 
report the best fit parameters and 90% confidence ranges of the model fit to the spectra in the three 
phases. We also list the result of fitting the flare spectrum subtracted from the pre-flare spectrum. 

The flare spectrum in 2009 shows a modest increase in temperature, changing from ~ 0.4 
keV prior to the event to ~ 0.9 keV (flare spectrum minus pre-flare spectrum). The temperature 
of the flare is similar to that observed in solar flares (~ 10 MK). The total flux of HD 189733A+b 
increases by a factor of 85% and the emission measure by 11%. The flux in the 0.3-8 keV band 
due to the flare is 1.3 x 10" 13 erg s" 1 cm -2 . Some features of the coronal spectrum are enhanced 
during the brightening. Fig. [5] shows the PN spectra in the three phases: black is pre-flare, red is 

o 

the flare and blue is after the flare. The most notable difference is around 0.57 keV or 22A where 
the O VII triplet is located (recombination: r = 21.62A, intercombination: i = 21.84A, forbidden: 
/ = 22.1 OA). This region of the spectrum remains in excess at a level of ~ 2<x with respect to the 
best fit thermal model. 

We explore this spectral region in details by means of the RGS 1 spectrum that allows a much 
higher resolution with respect to EPIC cameras in order to distinguish between the r, i, and / 
triplet lines of O VII. The relative intensity of these lines offers a diagnostic of the plasma density, 
wherein the forbidden line is formed when the plasma has low density and the corresponding 
energetic level is de-populated by collisions under increasing plasma density. Fig. [6] shows the 
spectra in the three phases as above. We have rebinned the spectra, by decreasing the instrumental 
resolution by a factor 5 to enhance the signal. In these spectra we can easily distinguish between r, 
i and / lines. During the flare the i line seems to disappear, and the / line is less luminous. After 
the flare the / line remains quite luminous and this fact suggests that low density plasma continues 
to be visible after the heating. There are indications that the r recombination line is shifted by 



> 15 km/s, which is consistent with the orbital speed of the planet (~150 km/s). Linsky et al. 



(|2010 ) looked for UV signatures of ionized gas co-orbiting with HD 209458b, while they found 
evidence of an exosphere, they concluded that higher S/N were needed to ascertain the velocity 
structure. Likewise, our data are not sufficient for such detection. Other visible changes during the 
flare which may be related to the increase of temperature are: an enhancement of O VIII line at 
18.96A Fe XVII lines at 15A and 17 A and an enhancement of the continuum level, 

In the framework of modeling of the flares in the stars in analogy with solar flares, several 
studies have been don e involving both simple models based on scaling laws and more detailed 



MHD simulations (see lRealdl2007L and references therein). A simple model for the flare is given 
by a loop formed by the magnetic field in which the plasma is constrained to move only along the 
field lines. The heating of the flare is provided at some point of the loop, which usually has an 
aspect ratio (base radius on loop semi-length) of r/L = 0.1 similar to loops observed on the Sun. 
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If the flare we observed took place in a coronal loop anchored to t he star surface, a r ough estimate 
of the semi-length of the loop is provided by the fo rmula given by ISerio et al.l (119 91). We obtain a 
semi-length of order of ~ 75% R* or 4 x 10 10 cm (|Serio et al.ll 1991c LReale et al.ll 19971) . assuming 
a peak temperature of 0.86 keV and no further heating mechanism at work after the initial loop 
heating. Given that we cannot estimate the true peak temperature, the size we would have derived 
is a lower limit. We can suppose that the size of the flaring loop is of the order of the stellar radius. 
Loops of the order of the stellar radius are inferred als o from a detailed modeling of a large flare 
observed in the very active dM star Proxima Centauri (IReale et al.ll2004|) . Given the proximity of 
the planet (at d ~ 8.75 R*) the perturbation on the stellar magnetic field due to the planet could 
affect the formation of loops and their evolution. The unusual flaring volume suggests that the site 
of the flare should be an extended region of the corona or an arcade formed of similar loops. From 
the estimates of the emission measure and the volume we derive a mean electron density in the 
plasma of n e ~ 9 x 10 9 cm -3 . 



From the ro tational period (P = 1 1.95<i. lHenry & Winnll2008l) and stellar radius (R = 0.767?*, 
Pont et al.l 120071) we infer a mean rotational velocity of v rot = 3.4 km s _1 . Early studies of MS 
solar type stars have s hown that the X-ray luminosity is well correlated with rotational velocity 
( Pal lavicini et al. 198 II), given the link between rotation, dynamo efficiency, coronal magnetic field, 
X-ray activity and L x . Comparing with the Sun we see that L x and v rot scale reasonably well with 
the L x ~ (vsinz) 2 law in the non saturated regime of L x . From the rotational velocity ratio between 
HD 189733A and the Sun (2 km s" 1 ) a mean magnetic field strength ~1.7 times the mean solar 
field is inferred. 

The M-type companion is not detected in either observations although it could have been 
resolved given the separation of ~ 12". We have estimated an upper limit to its X-ray luminosity 
in 0.3-8.0 keV band of L x < 9 x 10 26 erg s _1 . The non-detection of this star has some implications 
for the estimate of the age of the system as we will discuss in sect. |4j 



3.1. The X-ray Source near HD 189733 

The X-ray images in 2007 and 2009 show a source ~ 12" to the south of HD 189733 A. The 
source is brighter in 2007 than in 2009 (see Fig. [7]). The spectra obtained in 2007 are shown in Fig. 
[81 those from 2009 observation have poorer statistics and are not shown here. The source reveals a 
quite hard spectrum and, given that HD 189733A has a quite soft spectrum, the contamination from 
this source to HD 189733A spectrum is negligible. We have tried best fit with models composed 
by an absorbed power law and an absorbed APEC thermal model. In both cases the absorption is 
high: 7Y H ~ 1.4 X 10 21 cm" 2 . Given that the spectra from 2009 have poor statistics, a simple model 
with an absorbed power law of index of ~ 1 .7 has a good agreement with the data. On the other 
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hand, the spectra of 2007 have better statistics and the shape reveals some features around 1 keV 
so the simple power law is not appropriate. When modeled with an APEC thermal component 
we find a high plasma temperature value, of order of 8 keV. The flux in 0.3-8.0 keV band is 
1.9 x 10" 13 erg s" 1 cm" 2 in 2007, whereas it is 1.4 x 10~ 13 erg s" 1 cm -2 in 2009. The ne arest 
optical / in frared counterpart to t his X-ray source is an object in the 2MASS (ICutri et al. l2003h and 
NOMAD dZacharias et al]|2005|) catalogs at ~ 3". It has a poorly determined 2MASS photometry 
given the proximity to HD 189733 and quite large proper motions that hint that the cataloged 
object should be quite close. On the other hand, the value of absorption of the X-ray spectrum of 
this source suggests that it should be quite distant from the Sun but still inside the Galaxy. In fact, 
typical values for the total N H column in the direction of HD 189733 are between 3.8 x 10 21 cm -2 
and 3.9 x 10 21 cm" 2 (|Kalberla et al.ll2005l : iDickey & Lockmanl ll990). By considering the flux in 
2007 and 2009, and placing it at a distance of 1000 pc the luminosity yields logL x = 31.2-31.4. 



4. Discussion. 

In this paper we have analyzed two X-ray observations of HD 189733A: the first obtained 
during a planetary transit on April 17 2007 and the second during a planetary eclipse occurred on 
May 18 2009. Here we discuss the features observed in the light curve of 2009 and we relate them 
to possible effects of interaction of the planet and the host star. 



4.1. Observed features 



The first feature we have observed in 2009 is a softening of the emission during the planet 
eclipse at a level of > 3<x. Although similar variations of the spectrum hardness are seen in the light 
curve, this softening has a direct correspondence with the eclipse timing and, in addiction, during 
the eclipse the star has the softest emission even recorded in both 2007 and 2009 observations. The 
mean PI in EPIC-PN in the 0.3-1.5 keV band is 700 ± 10 eV before the eclipse and after the flare 
that occur post-eclipse, while it is 660 + 10 eV at the mid eclipse. 

An analogy between HD 189733b and Jupiter here may be appropriate. Jupiter shows en- 
hanced soft emission in X-rays coming from the equ atorial region and due to the inte raction of 
solar X-rays with the upper atmosphere of the planet (IBranduardi-Raymont et al.ll2007|) . Further- 
more, a harder X-ray emission is detected near the poles of Jupiter and may b e due to charge 
exchange mechanisms caused by precipitating ions of carbon, sulfur and oxygen (IHui et al.ll2009i : 
Branduardi-Ravmont et al. I I2007L The flux in 0.2-7 keV coming from Jupiter aurorae is ~ 7 x 10 
erg s" 1 cm" 2 that corresponds to a luminosity of ~ 5 x 10 15 erg s" 1 . By using simple scaling, taking 
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into account the star-planet distance in HD 189733, the higher luminosity of the star with respect 
to the Sun (10xL x , Q ) and the distance to HD 189733 we estimate that the power emitted by aurorae 
would be L Xi aur ~ 1.4 x 10 21 erg s" 1 , corresponding to a flux to Earth of fx,am ~ 3 x 10~ 20 , which 
is too faint to be detected by XMM-Newton, unless there are additional sources of energy caused 
by the magnetic field topology. If the softening is related to the secondary eclipse it cannot be 
related to a scaled emission of Jupiter for the case of HD 189733b. With only one observation we 
are not in the position to assess that the softening is related to the secondary eclipse. A series of 
observations at the same phase would reveal if the softening is systematically recurrent. 

The second relevant feature present in the light curve of 2009 is the most intense flare recorded 
in both observations, occurring ~3 ks after the secondary eclipse. The increase of flux is equal to 
~85% of the emission before the eclipse and the plasma temperature shows a modest heating 
passing from ~0.4 keV prior to the event to ~ 0.9 keV. The properties of the flare like temperature, 
luminosity, decay rate are similar to other flares observed in solar type stars without hot Jupiters. 
During the flare the lines of OVII at ~ 22A are more luminous than before and after the event. 
The OVII complex remains in excess with respect to the best fit model of the flare spectrum. High 
resolution RGS spectra show that the / OVII line remain quite luminous after the event and the i 
line almost disappears. 

During its cooling phase, we see several impulsive events (see Fig. [3]). It is as though there 
are multiple small flares occurring as the overall structure is cooling. The faster decay time would 
imply loops smaller than the main flaring loop. A peculiar characteristic of the small impulses 
is that the rate increases while the median PI steadily decreases. One way to understand this is 
to imagine that the flare occurred close to the limb, and at least in occurrence of such events the 
corresponding part of the corona was emerging into view. If our interpretation of this evidence is 
correct, then the flare is occurring on the limb when the planet is at phase approximately 0.54. This 
would put the flare displaced by ~ 77" and the meeting the sub planetary point (see figure HI). 

Given the timing of the flare one could speculate that the coronal region where it took place 
could be displaced by ~ 77" and trailing the sub-planetary point. The presence of spots and active 
regions on the stellar surf ace rotating synchrono usly with the planet have b een reported in literature 



in stars with hot Jupiters. iShkolnik et all (120031) and lShkolnik et all (120051) report that HD 179949 
and v And show variability of chromospheric activity indicators phased with planet rotation hinting 
thus the presence hot spots rotating synchronously with the planet as a consequence of star-planet 
interaction. The spots are displaced with respect to the sub-planetary point by +70° and -170" 
respectively. Analo gously to sol a r coro na these spots should be the chromospheric footpoints 



of coronal loops. IShkolnik et alj ((2008) has suggested also some excess of flaring activity on 
HD 189733 in phase with the planet rotation superimposed to the main long lasting active regions 
on the stellar surface. 
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One speculative hypothesis is that the large flare we have observed could arise from active 
regions for which the activity is triggered by the complex magnetic interaction that should be 
present between the star and the planet. In summary, we cannot draw a conclusion on this point 
with only one observation. As in the case of the softening, a series of further observations around 
the secondary eclipse could establish the systematic recurrence of high activity just after the end 
of the eclipse and coming from active regions triggered by star-planet interaction. 



4.2. Comparison with SPI models 



Lanzal (|2008|) proposed an analytical model for the stellar magnetic field in the case of SPI. 
This model was aimed to explain the phase shifts of active regions with respect to the sub-planetary 
point observed. With a suitable choice of the model parameters the author shows that this simple 
model can reproduce these phase shifts for HD 189733 HD 179949, v And and r Bootis. The 
prediction from this model for the case of HD 189733 is that a p hase sh i ft betw een sub-planetary 
point and active region should be ~ 78°. With the same model Lanza ( 2009 ) has addressed the 
energy budget, the intermittent nature of the visibility of spots and the effect that could be seen in 
radio and X-ray bands. From a qualitative point of view the stellar magnetic field should be able 
to drag material evaporating from the planet and to deposit it in a azimuthal flux rope in the outer 
part of the corona. This region should be the site of an enhanced X-ray activity due to magnetic 
reconnection (see Fig. @]for a simple schematic). By comparing our observational results to this 
model, we observe that the la rge flare in 200 9 can fi t consistently with this model and the excess of 
flaring activity suggested by iShkolnik et al.1 (|2008|) under the hypothesis that it come from active 
regions displaced by ~ 75° - 78° with respect to the sub-planetary point. 

We have carried out prel iminary MHD simu lations of the star+planet system within the dipo- 
lar field setup ("Case A") of Cohen et al.l (|2009|) . but with the stellar and planetary parameters of 
the HD 189733 system. The planetary magnetic field is taken to be half that of Jupiter (with an 
equatorial field of 2 G); this is a reasonable assumption since tidal effects will have reduced the 
planetary rotation compared to that of Jupiter. Figure [9] shows one frame of the time evolution of 
the magnetic field of HD 1 89733A+b system based on our MHD simulation. The simplified dipole 
model produces a density enhancement in the intervening interplanetary medium that is threaded 
by complex magnetic fields. There is thus a greater potential to observe transients occurring in this 
region. The density enhancement is formed by plasma trapped into closed magnetic field lines that 
would have been opened in the absence of the planet interaction. 

In summary, we find elements of consistency between the observed features in X-ray band, 
optical chromospheric indicators and model predictions. The magnetic interaction of planet and 
host star can enhance the activity in some regions of the stellar surface and that probably these are 
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among the most active regions of the star chromosphere and corona. As suggested by lShkolnik et al 



(|2008|) basing on chromospheric activity indicators, the spot activity can be enhanced or reduced 
due to the complex magnetic field evolution. This part of the star-planet interaction is observed as 
flares. On the day side of the planet, it is plausible that the interaction manifests itself as colliding 
wind from the star on to the inflated atmosphere and with polar aurorae that emit X-rays harder 
than the main coronal emission, although this emission remains too faint to be detected. 

Although our data do not allow us to draw specific conclusions, they suggest that effects of 
the magnetic interaction between a star and a close planet could be detected in X-rays. A series 
of observations in the same phase range near the secondary eclipse and perhaps coordinated with 
optical observations of chromospheric activity indicators can allow to distinguish systematic effects 
caused by planet presence from random variability intrinsic to the coronal activity. 



4.3. The age of HD 189733 and its activity 



The age of HD 189733 has been estimated bv iMelo et aP foood) to be greater than 0.6 Gyr 



basing on indicators of chromospheric activity calibrated on normal Main Sequence stars. Fur- 
thermore these authors reported a very low abundance of Li (A(Li)< -0.1) in HD 189733A. For 
comparison in the Sun A(Li) is ~ 2 at an age of 4.6 Gyr. Lithium is destroyed during Main Se- 
quence phase because of convection that brings lithium in the inner layers of the star where at a 
temperature of 2-3 MK is burned. Both Li and chromospheric indicators are sensitive for ages 
< 0.5 - 1 Gyr. In the case of HD 189733 the use of chromospheric lines could be biased by the 
enhanced activity of the star itself and, eventually, by SPI. Our X-ray observation suggests a better 
constraint on its age. 

The non-detection of HD 189733B, a M4 spectral type star, at a level of L x < 10 27 erg s _1 
suggests an age estimate of HD 189733 older than the current one. M-type stars exhibit a decay of 
X-ray emission during Main Sequence phase slower than in G and K type stars due to the lower 
rate of dissipat i on of angular momentum and the link between rotation and dynamo efficiency. 
Feigelson et al. d2004h have studied a sample of low mass stars in the Chandra Deep Field North. 



In this sample M4 stars are characterized by L x ~ 10 28 erg s _1 , and the age of the sample is 
presumed to be greater than 3 Gyr. Given the non-detection of HD 189733B we suggest that the 
age of this system could be tentatively estimated older than 2 Gyr. 

Placing the system at an age of 2 Gyr the X-ray luminosity of HD 189733 (L x ,« ~ 10L x ,s U n) 
fits the trend observed from studies of X-ray l uminosities and activity in K-type stars of open 



cluste rs and field stars (see iGiardino et al.ll2008l and references therein). iMamajek & Hillenbrand 



(|2008|) consider the age estimates from various indicators like Ca H&K lines, A(Li), rotation and 
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X-ray emission, obtaining thus several relations between age and these observables. With a ro 
tational period of ~ 1 2d and a Kl type star, HD 189733 would fit close to the age of Hyades 



(Ma majek & Hillenbrandll2008l see their fig. 10). But the value of logR' HK reported by lMelo et al 



(I2006|) and the curves given by iMamajek & Hillenbrandl (I2008|) in their fig. 12 would suggest an 
age of r > 1.7 Gyr. Furthermore, U, V, W space motion components (U = -8.6 km /s, V = -14.4 



km/s ? W = 015.9 km/s) are not compatible with those of Hyades moving group (|Montes et al. 



200 1|) . In conclusion, although the determination of the age of an isolated star can be difficult and 
poorly constrained, the non-detection of the M type companion below log L x ~ 27 hints an age of 
HD 189733 around 1.7-2 Gyr. 



5. Conclusions 

In this paper, we have presented the analysis and the interpretation of X-rays light curves and 
spectra that we have acquired during a secondary eclipse of the planet HD 189733b observed in 
2009 with the XMM-Newton satellite. For completeness we have also analyzed in a homogeneous 
way a previous observation taken in 2007 during a planetary transit on the same target. Our group 
has carried also MHD simulations of the magnetic interaction between star and planet. From the 
observations we have estimated that the average X-ray flux of the star in 2009 (fx ~ 3.7 x 10~ 13 
erg s _1 cm" 2 ) has been has been 45% higher than the average flux in 2007. The light curve of 2007 
shows impulsive variability which is considered quite usual for a K-type star quoted at an age of 
1 Gyr. During the planetary transit we do not observe spectral or rate variations timed with the 
transit. 

In the light curve of 2009 observation we have observed a softening of the spectrum significant 
at a 3<t level during the secondary eclipse and the biggest flare recorded in both observations 
occurring 3 ks after the end of the planet eclipse. The softening is well phased with the planet 
eclipse. The median of the spectrum change from an average value of 700 + 10 eV before and 
after the eclipse to 650 + 10 eV at the eclipse center. The decay phase of the main flare in the light 
curve shows several superimposed minor ignitions perhaps occurring in the same region where the 
main flare took place. We have fitted the spectrum subtracted from the pre-flare spectrum with a 
mono-thermal APEC model. The flare exhibits a average plasma temperature of 0.9 keV while 
the quiescent plasma temperature is around 0.5 keV. The RGS spectrum shows a variability of the 
relative strength of the lines of OVII triplet prior, during and after this flare. In particular we noticed 
the reduction of inter-combination and forbidden lines during the flare, and a bright forbidden line 
after the end of the flare. These two behaviors of the OVII triplet lines are explained as an increase 
of density in the flaring region during the ignition, and a persistent low de nsity bright/hot pl asma 



after the end of the flare. By using a simple model of cooling loop as in ISerio et alj (| 199 11) and 
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( 199 7J), we infer a size o f the flaring region of order of the stellar radius. From L x 



Reale et al. 

scaling law (|Pallavicini et al.lll981|) we infer a magnetic field twice the strength of the solar one. 



In order to explore the possibility that the planetary magnetosphere can directly affect the 
stellar coronal structure, we have conducted MH D simulations of a model of interacting magnetic 
fields with opposite alignment as in Case A of Icohen et al. ( 2009 ). From these simulations we 
infer that due to the relative motion of the planet and the star the global magnetic field is distorted 
and has an enhanced probability of reconnection events that can occur on the space between planet 
and star. Given the coupling of magnetic field and plasma, an enhanced plasma density is present 
around the planet and it forms a trail rotating synchronously with the planet. Globally we find 
consistency between observations and predictions from MHD simulations, hinting that the overall 
activity could be enhanced by magnetic SPI. Hence, the age estimates of the system using activity 
indicators are suspect. Currently the age of HD 189733 is estimated to be at least > 0.6 Gyr. The 
M-type companion, HD 189733A, is not detected at a level of luminosity of L x < 10 27 erg s _1 
suggesting thus an age for the system older than 1.7-2 Gyr. 

We have also analyzed the spectrum of a X-ray source close to HD 189733 and brighter 
in 2007 than in 2009. The spectrum is quite hard and from best fit we derive an absorption of 
/Vh ~ 1.4 x 10 21 cm -2 suggesting that this is a object in the Galaxy but quite distant from the Sun. 
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Fig. 7.— EPIC and 2MASS images of the X-ray source at ~ 12" from HD 189733. From left 
to right: MOSl > IkeV in 2007, MOS 1 > 1 keV in 2009, 2MASS J-band image. Positions of 
HD 189733A, HD 189733B and the source are indicated by crosses. 



Table 1: Best fit parameters of spectra collected during the observations in 2007 and 2009. For 
2009 we list the results of the three time intervals: before flare, during flare and after the flare. 
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Fig. 8. — PN and MOS spectra of the X-ray source at ~ 12" from HD 189733obtained during the 
2007 observation. The best fit model is an absorbed APEC thermal model. 
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Fig. 9. — MHD model of the star+planet. On the left panel: one frame of the MHD simulation. 
The star and the planet are represented by the red and blue spherical shades, respectively. Lines 
of the magnetic field are drawn in gray, the inset plot in the right bottom part show the central star 
and the planet. On the right panel: density image viewed from the top. Model predict a denser trail 
orbiting synchronously with the planet. 



